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GLAUBER MODEL FOR HEAVY ION COLLISIONS FROM LOW 
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The Glauber model is extensively applied to heavy ion collision for describing a num¬ 
ber of interaction processes over a wide range of energies from near the Coulomb barrier 
to higher energies. The model gives the nucleus-nucleus interaction in terms of interac¬ 
tion between the constituent nucleons with a given density distribution. The model is a 
semiclassical model picturing the nuclear collision in the impact parameter representation 
where the nuclei move along the collision direction in a straight path. In these lectures 
we derive this model and discuss its applications in variety of problems in nuclear and 
high energy physics. 

I. CROSS SECTIONS AND IMPACT PARAMETER REPRESENTATION 

The total reaction cross section in the collision of two nuclei as per the partial wave 
analysis is given by 

00 

'yR = TsllC^i + m-\Si?)- ( 1 ) 

^ i=0 

Here Si = exp(—2*5/) is called the scattering matrix and 5i is the nuclear phase shift. 
The factor (1 — |S';p) is called the transmission coefficient and |S';p is referred to as 
transparency function or the probability that the projectile undergoes no interaction at a 
given /. In a semiclassical approximation, one can write angular momentum I in terms of 
momentum k and impact parameter b as 
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/ + i = kb. (2) 

Thus, using (2/ + 1) = 2kb andX); = k J db we can get the reaction cross section in impact 


parameter representation as 

aR = 27r J bdb{l-\S{b)\‘^). (3) 

The factor Tr{b) = 1 — |5'(6)p is nothing but the transmission coeff. If the two nuclei are 
assumed to be sharp spheres with radii -Ri and R 2 then 

Tr{b) = 1 for b < Ri + R 2 

= 0 for b > Ri + R 2 . (4) 

The total reaction cross section in this case is given by 

PR1+R2 

aR = 2 ti / bdb = 7i{Ri + R 2 )^. (5) 

Jo 

This is the well-known geometric formula for the cross section. 

The nucleus-nucleus differential cross section as a function of center of mass scattering 
angle 9 is described as 

where the scattering amplitude f{6) for the non identical spinless nuclei is given by partial 
wave analysis which can be written as a sum of Coulomb and nuclear scattering amplitudes 
as 


f{0) = fc{e) + fN{0) 


= fc{0) + ^ E(2/ + - l)Pz(cos0). 


(7) 


Here fc{9), the Coulomb scattering amplitude and ai, the Coulomb phase shift are given 
by the following expressions: 

O' 


f (a\ 0 20 

fc{0) = 2 


2 iao — 2 iri In sin 


77 

= aiiri) +tan~^{ 


i + r 


( 8 ) 


(9) 


The S wave Coulomb phase shift uq can be set equal to zero without any loss of generality. 


2 



II. THE GLAUBER MODEL 


The Glauber model of multiple collision processes provides a quantitative consideration 
of the geometrical conhguration of the nuclei when they collide. The Glauber model 
basically describe the nucleus-nucleus interaction in terms of elementary nucleon-nucleon 
interaction. It is based on the assumption that the nucleus travels in a straight line 
path. At high energies this approximation is very good. At low energies the nucleus 
is deflected from straight line path due to Goulomb repulsion. The so called Goulomb 
modified Glauber model (discussed in the next section) also gives very good description 
of heavy ion scattering at low energies. 

Gonsider the collision of a projectile nucleus A on a target nucleus B. Dehne t{h)dh as 
the probability for having a nucleon-nucleon collision within the transverse area element 
dh when one nucleon is situated at an impact parameter b relative to another nucleon 
which is normalized according to 

J t(h)dh = 1. (10) 

We dehne the probability of hnding a nucleon in the volume element dh^dzA in the nucleus 
A at the position (hA,ZA) is Pa(Joa, ZA)dhAdzA which is normalized as 

J pA{hA, ZA)dhAdZA = 1 . ( 11 ) 

Similarly, the probability of hnding a nucleon in the volume element dhsdzB in the nucleus 
B at the position {hBiZs) is pb^^B) ZB)dhBdzB which is normalized as 

[ Pb(^b, ZB)dhBdzB = 1 . ( 12 ) 
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FIG. 1. Collision of two nuclei at an impact parameter b 


The probability for occurrence of a nucleon-nucleon collision [see Fig. (|l])] when the 
nuclei A and B are situated at an impact parameter b relative to each other is given by 

T{b)aNN = j Pa{^a, ZA)dhAdzA PbO^b, ZB)dhBdzB t{h - bA + be) a^N- (13) 

This can be written in terms of z-integrated densities as 

T{b)aNN = J pAi^A)dhAp%{hB)dhBt{h - Ta + bg) aNN- (14) 


Here a^N is the inelastic nucleon nucleon cross section. Thus, the collision probability 
we are talking about is for an inelastic collision. There can be upto Ax B collision. The 
probability of occurrence of n collisions will be 


F(n, b) = 


'AB^ 


n 


(1 


(^) 


AB-r 


(15) 


Here, s = 1 — T{b)aisfN. The total probability for the occurrence of an inelastic event in 
the collision of A and B at an impact parameter b is 


™ n=l 


-,AB 


(16) 


4 











The total inelastic cross section can be written as 


<^.1® = 2,r IM6 (l - . (17) 

From here one can read the scattering matrix as 

\S{b)\^ = s^^ = {l-T{b)a^N)^^. (18) 

In the optical limit, where a nucleon of projectile undergoes only one collision in the target 
nucleus 

|5(6)|^ ~ exp{-T{b)aNNAB). (19) 

The scattering matrix can be dehned in terms of eikonal phase shift x(6) as 

S{b) = exp{-ix{b)). (20) 

Comparing Eq. ([T9| ) with Eq. (pO]), the imaginary part of eikonal phase shift is given by 

Imxib) =T{b)aNNAB/2. ( 21 ) 

If the ratio of real to imaginary part of NN scattering amplitude is ajsfN then real part of 
X{b) is 

Rex(6) = T{b)aNNO'NN^B/2. (22) 

Once we know the phase shift and thus the scattering matrix, we can calculate the reaction 
cross section and the angular distribution. 

III. CALCULATION OF T{B) IN MOMENTUM SPACE 

In the co-ordinate space T(6) is derived as 

Tib) = I pji(b^)db^ p%ihB)dhB t(b - bA + bB). (23) 
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It is a four dimensional integration: two over bA and two over be- It is convenient to 
write it in momentum space as 


T{h) = — / p\{hA)dhA PB{hB)dhB exp (-iq.(b - bA + be)) fNN{q)d^q- 


Here fwNiq) is the q dependence of NN scattering amplitude given by 


t(b) = 


/ e-^^-^fNNiq)d^q 


T{h) = J exp(-iq.b) ^^(bA) exp(iq.bA)dbA p^bs) exp(-iq.bB)dbij fNN{q)d^q 

= J e-^^-'^SA{q)SB{-q)fNN{q)d‘^q 

= ^ / Mqb)SA{(l)SB{-(l)fNN{q)qdq. (26) 

Here SA{q) and SB{—q) are the fourier transforms of the nuclear densities. The function 
fNN{q) is the fourier transform of the prohle function t{h). The prohle function for the 
NN scattering can be taken as delta function if the nucleons are point particles. In general 
it is taken as a gaussian function of width ro as 


t(b) = 


exp(— 


Thus, 


/iV7v(g) = J e"'^-^t{h)dh 

= ^ / e''’-^exp(-6Vr2)db 
Trrg J 

= expl-r^q^/A). (28) 

Here, rg = 0.439 fm^ is the range parameter and ctaat is the nucleon-nucleon inelastic 
cross section which is taken as 3.2 fm^ at high energies. 

A. Calculation of T(6) using Gaussian densities 


If the nuclear densities are assumed to be of the Gaussian shape given by 
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Pii.r) = pi(0)exp(-r /aj 


(* = 1 , 2 ). 


(29) 


Here the parameters Pi(0) and a* are adjusted to reproduce the experimentally determined 
surface texture of the nucleus. The 2 :-integrated density will be 


Piir) = pi{0)V^aiexp{-b‘^/a^i). 


This one can write in the momentum representation as 


Si{q) = Pi{0){^/waif exp{-q^ay A). 


The overlap integral T{b) can be written as 

1 


T(6) = 


(27r)2 

1 

(27r)2 


e *‘^-'’5'i(q)S'2(-q)/iVAr(g)d^g 


7r^pi(0)p2(0)aia2 J e exp(—a^g^/4)(i^g, 


where + Tq. Performing q integration we get 


^(^) = 7T^^Vi(0)p2(0)a?a^ 


= TT 


(27r)2 

2Pi(0)p2(0)afa| 


Att 


a 

,2 // 2 


exp(-6^/a ) 


of + ^2 T ^0 


exp(-6 /(a^ + 02 + ro)). 


(30) 


(31) 


(32) 


(33) 


B. Calculation of T{b) using 2pf densities 


The two parameter fermi density is given by 


p(r) = 


do 


l + exp(^)’ 

where po = 3/ (47rc^{l + ^^^))- Thus, the momentum density can be derived as 

Sttpo f . z{l + e-2^) ^ 

Siq) = ^3 i_e- 2 . (^sma; ^ - xcosa ;j . 

Where = ndq and x = cq. Here d is the diffuseness and c, the half value radius in 
terms of rms radius R for the 2pf distribution is calculated by c = (5/3i?^ — 7. 
The equation can be solved numerically for this density and the overlap integral can be 
extracted. 


(34) 


(35) 
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C. NN scattering parameters cttvat and 


In this section, we give the expressions for and obtained from the parame¬ 
terization of experimentally measured cross section. The average aisf^ in terms of proton 
numbers {Zp and Zp) and neutron number {Np and Np) of projectile and target nuclei 
is written as 


O'NN — 


NpNpann + ZpZpapp + (ZpNp + NpZp)a, 


np 


ApA. 


p/ip 


(36) 


Here, pp cross section app and nn cross section ann are given in (fm^) by 


app = (Tnn = 1-373 - 1.504//3 + O.STG//^^ + 6.867/31 


(37) 


Here (3 is the velocity of projectile nucleon. For np cross section two expression are 
used. If the energy per nucleon Ejq > 10 MeV, then 


anp = -7.067 - 1.818//3 + 2A2Q/(3‘^ + 11.385/3. 


For Epi < 10 MeV, 


273 


1763 


^np 


(1. - 0.0553.Fjv)2 + 0.35.FJV (1- + OASAEn)^ + 6.8En ' 

The average is written as 

NpNpOirin^nn T ZpZpUppCpp (^ZpNp -|- NpZp^CXnp^^np 

= a^^ApAp • 

The parameterized forms of app, ann and Unp are written as 

app = ann = 0.1810 -|- 4.0818p -|- 0.3327p^ 


and 


ttnp = —.698 -|- 4.9762p — 1.277p^, 


(38) 


(39) 


(40) 


(41) 


( 42 ) 


where p is the momentum of projectile nucleon in GeV. 



IV. COULOMB MODIFIED GLAUBER MODEL 


In the nucleus-nucleus collision at low energies it is the surface region of the nucleus 
that contributes to the scattering amplitude in a non-trivial manner. Since V/E « 1 
even for low and intermediate energies in the surface region, the eikonal approximation 
can be extended to the low energy heavy ion collisions. The basic assumption of the 
Glauber model is the description of the relative motion of the two nuclei in terms of 
straight line trajectory. The density overlaps are evaluated along straight lines associated 
with each impact parameter b. The modihcation in the straight line trajectory due to 
the Coulomb held can not be ignored, especially in the case of heavily charged systems 
at relatively low bombarding energies. 
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FIG. 2. The straight line trajectory is assumed at the distance of closest approach 


For low energy heavy ion reactions the straight line trajectory is assumed at the 
distance of closest approach Vc [see Fig. (^] calculated under the influence of the Coulomb 
potentials for each impact parameter b as given by, 

Tc = [r] + (43) 


which is a solution of the equation 


E- 


r 


2fi r2 


0 , 


62 

^ 2 


0 , 


kr^ — 2r/r — b^k — 0. 


(44) 

(45) 

(46) 


Here )/ is the dimensionless Sommerfeld parameter defined as 

Z\Z^^ k 

^ ^ hv ^ E 2' 

Thus hnally the reaction cross section will be 


(47) 


AB _ 
^in ~ 


y M6 (l - (1 - T{rc)aNK)'*^) ■ 


(48) 


Figure (^) shows the reaction cross section for O + 2 °^ Pb system as a function of 
projectile energy per nucleon calculated with Glauber model (GM) and Coulomb modihed 
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Glauber model (CMGM) along with the data. At higher energies both GM and GMGM 
result merge. But at lower energies one requires to take Goulomb modihcation in the 
trajectory. 



FIG. 3. Reaction cross section for O + Pb system calculated with Glauber model 
(GM) and Goulomb modified Glauber model (CMGM) along with the data 

V. PARTICIPANT-SPECTATOR PICTURE IN HIGH ENERGY HEAVY ION 

COLLISIONS 

Nucleus-nucleus collisions at ultrarelativistic energies are being explored in order to 
search for the formation of quark gluon plasma. Already several experimental results 
are available. One of the stategies is to describe these results as the superposition of 
the nucleon-nucleon collisions and look for possible departures from this prescription. 
When two nuclei collide at these energies, the nucleons which come in the overlap region 
depending on the impact parameter are called participant nucleons and those which do 
not participate are called spectators, [see Fig. (I)] 


11 






Spectators 



These spectators and participants nncleon decide how mnch energy is going in the 
forward and how much in the transverse direction which can be measured. From these 
measured energies one can know the impact parameter for a particular event. Let us write 
the total inelastic cross section for a nucleus-nucleus collision as 

= ^TT y bdb (l - (1 - T{b)aNN)'^^) ■ (49) 

The term (1 — T{b)a]^]\f)^^ gives the probability that in a nucleus nucleus collision none 
of the nucleons collided with each other. For nucleon-nucleon collision A = B = 1 thus 
s = (1 — T{b)aNN) gives the probability that two nucleon at an impact parameter b do 
not collide. 

When two nuclei A and B collide the probability of nucleon remaining in the projectile 
will be 

Pp = s{bf (50) 

and the probability of nucleon remaining in the target will be 

Pt = s{b)^. (51) 
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The probability of having a participant nucleons from the nucleus A is given by bino¬ 
mial distribution as 

P(a,l,)^A'j(l-s‘>r(sT-‘- (52) 

Similarly, 

m5)= (53) 

The number of collisions will be 

P(n,6)= i'^®V-s)”(i,)“-” (54) 

The average number of projectile participant and its standard deviation for each impact 
parameter b is given by 

< a >= A[1 — s{b)^], (55) 

al = A[1 - s{b)^]s{b)^. (56) 

The average number of target participant and its standard deviation for each impact 
parameter b is given by 

< (3 >= B[1 - s{b)% (57) 

al = B[1 - s{b)^]s{b)^. (58) 

The average number of total participant from both projectile and target for each impact 
parameter is given by 

-^participants = A[l — s{b)^] + B[l — s(5)"^]. (59) 

The average number of N-N collisions is given by 

Ncollisions = AB[1 - s(6)]. (60) 
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The Fig. (^) shows the number of participants and the number of N-N collisions as 
a function of impact parameter with the expressions given in the text along with the 
calculations with the Monte Carlo code FRITIOF and Geometric model of J. Cosset et 

ah 




FIG. 5. (a) The number of binary collisions as a function of impact parameter calculated with 
the expressions given in the text along with the Monte Carlo code FRITIOF and Geometric model 
of J. Cosset et al. (b) The total number of participants as function of impact parameter. 


A. Forward energy 


The forward energy is measured in the detector called zero degree calorimeter (ZDC). 
The energy measured in the ZDC is can be related to the number of spectators as 


Ez{oi) — {A — a)Eo, 

where Eq is the projectile energy per nucleon in the lab frame. 
The cross section of a nucleons participating from projectile is 


aa = 27r J P{a, b)bdb. 


The forward energy flow cross section then can be written as 

da (To 
dEz Eq ’ 


(61) 


(62) 


( 63 ) 
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taking into account the discrete nature of the variable a. The average forward energy for 
a collision is written as 


Ez{a) = Eo{A- <a>) = EoA[l - (1 - s'®)] = EqAs^. (64) 


B. Excitation energy 

When the nucleons collide at relativistic energies they become excited nucleons which 
due to time dilation become long-lived and decay into secondary fragments outside the 
nuclei. This energy which goes in exciting the nucleons is called excitation energy which 
manifests in the measurement of transverse energy and multiplicity. The (maximum) 
excitation energy which is related to the number of participants is written as 

Egx Efyui -f- /3). (65) 

The centre of mass energy for this case is given as 

El^ = {{aEA + PEsf - {aPA + PPs?), ( 66 ) 

where 

Pa = \lE\-m% and Pb = \jE^ - m%. (67) 

The cross section for having a participants from A and (3 participants from B is 

aaf3 = 2% j P{A,a,b)P{B, (3,b)bdb. (68) 

Here Ea and are the total energy and momentum of the nucleus A per nucleon and 
mjsf is nucleon mass. The cross sections for all the possible combinations of a and (3 are 
calculated and the results are binned to obtain the cross section for each of the excitation 
energy bin. 
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C. Multiplicity 


The total charged particle multiplicity in pp collision for a center of mass energy ^/s 
is given by 


< ua >= 0.88 + 0.44 Iri(s) + 0,118(lri(s))^ 


( 69 ) 


The center of mass energy for pair of participants from nucleus A and nucleus B can be 
written as 


^af3 



(70) 


Thus, the total multiplicity in nucleus-nucleus collision in the case of a nucleons partici¬ 
pating from nucleus A and (3 participating from nucleus B can be written as 

<M>= (0.88 + 0.44 + 0A18{\n{sa^)f) ■ (71) 

The cross sections and the multiplicity for all the possible combinations of projectile and 
target participants are calculated and the results are binned to obtain the cross section 
for each of the multiplicity bin. 


VI. J/V' PRODUCTION CROSS SECTION 


The probability for j/'ip production in a nucleus nucleus collision can be found out if 
we know the j/p) production cross section in a nucleon nucleon collision, 

= ABafllnb), (72) 

The cross section is 

I dbT(t,). (73) 

The integral is 1 as per the normalization condition 
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(74) 


AB 


a 


NN 

i/V" 


= AB. 


The j/ijj once produced can collide with the other nucleons present in the nucleus and 
thus it is possible that it is absorbed (converted into D{cu) and D{cu)) during the path 
it travels in the nucleus. This nuclear absorption can be taken into account by writing 


a. 


AB 


= AB exp{-Laabspo)- 


NN 


(75) 


Here, po = 0.14 fm ^ is the equilibrium nuclear matter density, aabs = 6.2 mb is the 
absorption cross section and L is the average path length given by 


L — La + Lb — - ^— Ra H-^— Rb^ , (76) 

where, Ra = r^A^R and Rb = tqB^R. 

Figure (P) shows the suppression of J/T over Drell-Yan pairs as a function of the 
average nuclear path length L of the cc pre-resonance, for NA38, NA50 and NA51 data. 



FIG. 6. The suppression of J/'h over Drell-Yan pairs in 2.9 GeV < M < 4.5 GeV as a 
function of the average nuclear path length L of the cc pre-resonance, for NA38, NASO and 


NA51 data. The NASO Pb-|-Pb data (filled circles) is from the 1996 run. 


17 








VII. THE MONTE CARLO EVENT GENERATORS 


There are many Monte Carlo modes available to study the relativistic heavy ion colli¬ 
sions as superposition of nucleon-nucleon collisions such as FRITIOF and IRIS. In these 
codes the two nuclei are made to interact at a random impact parameter. A Woods-Saxon 
distributed nuclear density, giving the probability of hnding a nucleon at a given distance 
from the centre of the nucleus, is used. Each nucleon is given space coordinates and 
the number of subcollisions is calculated letting a projectile hit the nucleus at a random 
impact parameter. Here a frozen straight line geometry is used and all nucleons inside a 
cylinder surrounding the path of the projectile are considered to participate. The num¬ 
ber of binary collisions is recorded and the two partners involved in each subcollision are 
determined. In each subcollision, momentum is exchanged. 

FRITIOF is based on the LUND string picture of hadron-hadron interactions. In this 
picture, each nucleon-nucleon collision results in excitation of the nucleon by the stretch¬ 
ing of a string between the valence quark and diquark. In nucleus-nucleus interactions 
each nucleon can make several encounters, the objects can get further excited, thereby 
increasing their masses during the passage through the nucleus. A phenomenological ex¬ 
citation function determines the mass and momentum of the string after each interaction. 
Finally, all the excited objects hadronize independently, like massless relativistic strings, 
according to the Lund model for the jet fragmentation. The hadronization takes place 
outside the nuclei and thus no intranuclear cascading is considered. 

Three other models based on string picture of hadron-hadron interactions are IRIS, 
MCFM and VENUS all of which are colour exchange models based on the Dual Parton 
Model (DPM) of Capella et al. Here the basic mechanism of string formation is colour 
exchange between the quarks of the colliding nucleons. In these models the string prop¬ 
erties can be calculated from structure functions. HIJET is an extension of the ISAJET 
model of hadron interactions and MARCO is based upon a phenomenological parameter¬ 
ization of nucleon-nucleon collisions. ATTILA has the possibility of rope formation from 
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overlapping strings and MCFM and HIJET which allow cascading of produced particles 
when the assumed particle formation time is short. RQMD is an extension to relativistic 
energies of the Quantum Molecular Dynamics which has been applied to nucleus-nucleus 
collisions at much lower energies. 

The nucleus-nucleus collision geometry should be the same for all of the Monte-Carlo 
codes since nuclear density distributions are well-known from nuclear physics and should 
not be treated as free parameters. It should be noted that deviations as small as 5-10 % in 
the treatment of the nuclear geometry are signihcant since one hopes to draw conclusions 
about deviations from the measurements which are of similar magnitude. It is clear 
that before hrm conclusions can be drawn upon the signihcance of the differences in the 
physics of the models at the nucleon-nucleon level, or on whether there exists experimental 
evidence for rescattering of the produced particles, or even for QGP formation, it will be 
necessary to ensure that the various models treat the nuclear geometry correctly and 
consistently. 
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